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4^ Abstract 
o 

g We derive the criterion for the Bose-Einstein condensation (BEC) of a Gaussian field (p (real or 

^ complex) in the thermodynamic limit. The field is characterized by its covariance function and 

the control parameter is the intensity u = Hv?] 12/^1 where V is the volume of the box containing 
the field. We show that for any dimension d (including d = 1), there is a class of covariance 

g 

I functions for which exhibits a BEC as u is increased through a critical value Uc- In this case, we 

C investigate the probability distribution of the part of u contained in the condensate. We show that 

I ^1 depending on the parameters characterizing the covariance function and the dimension d, there can 

04 be two distinct types of condensate: a Gaussian distributed "normal" condensate with fluctuations 

> ^ 

0\ scaling as 1/vF, and a non Gaussian distributed "anomalous" condensate. A detailed analysis 

^ of the anomalous condensate is performed for a one-dimensional system [d = 1). Extending 

^ this one-dimensional analysis to exactly the point of transition between normal and anomalous 

condensations, we find that the condensate at the transition point is still Gaussian distributed 



^ but with anomalously large fiuctuations scaling as ■\/\n{L) / L, where L is the system length. The 

^ conditional spectral density of knowing u, is given for all the regimes (with and without BEC). 
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I. INTRODUCTION 



This paper is devoted to the concentration properties of a Gaussian random field </? (real 
or complex) in the "thermodynamic" limit, V — >■ +00 with fixed intensity u — Hv^HI/^, 
where V is the volume of the box containing the field and \\(p\\l = f^^^ \(p{x)\'^ d'^x. For a 
finite system (with given V < +00), the question was first addressed in [1] heuristically 
and numerically in the context of laser-plasma interaction physics with spatially smoothed 
laser beams and, more generally, of linear amplification in systems driven by the square of 
a Gaussian noise. The results of [1] were extended and given a mathematically rigorous 
meaning in [2]. 

Write C^^* the covariance of </? in a finite box with given V < +00, and ki the largest 
eigenvalue of C^^*. The main result of [1, 2] states that (p concentrates onto the cigcnspace 
associated with ki as \\'p\\2 gets large. By revealing that the realizations of a Gaussian 
field that may cause the breakdown of a linear amplifier are delocalized modes, this result 
overturned conventional wisdom ^ that breakdown is due to localized high-intensity peaks of 
the Gaussian field (the so-called high-intensity speckles, or hot spots [3]). 

One of the goals of this paper is to study and provide a physical interpretation of this 
concentration property of a Gaussian field in the thermodynamic limit. We show that the 
emergence of a delocalized mode as the intensity u increases beyond a critical value Uc is 
similar to the Bose-Einstein condensation in an ideal Bose gas: the eigenspace associated 
with the largest eigenvalue Ki plays the role of the ground state in the Bose gas. Technically, 
this is also very similar to the spherical model of a ferromagnet. In that model, the nearest- 
neighbor ferromagnetic interaction of continuous spins Sj e R on a lattice is supplemented 
with the global constraint Y^i-^i = ^> where N is the number of lattice sites [4]. For 
large N and in three or higher dimensions, the Fourier component of the spin field with 
A; = gets thermodynamically populated as one reduces the temperature below a critical 
temperature T^., signalling the onset of a global nonzero magnetization (a delocalized mode). 
The mechanism of this phase transition is similar to the Bose-Einstein condensation [5] . In 
our problem, the role of the constraint sf — N is played by \\(p\\l — J^^^^ Ivi^)]"^ d'^'x — uV , 

^ Most theoretical models dealing with stochastic amplification beyond the perturbative regime are hot spot 
models. In view of this result, the implicit assumption about the leading role of hot spots underlying all 
these models should be carefully reexamined. 
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with fixed u. Our work thus makes a nice hnk between the concentration properties of 
random fields with fixed intensity (of much interest in e.g. laser-plasma physics) and the 
traditional Bose-Einstein condensation in statistical and atomic physics. Indeed, it can be 

shown that for u large enough, the only contribution of the eigenspace associated with k,i is 
greater than the one of all the other eigenmodes of C^^* (which remains bounded) [2] . Such 
a behavior is typical of a Bose-Einstein condensation onto the eigenspace associated with 

Of course, like any other phase transition, no sharp condensation can occur in a finite size 
system. The thermodynamic limit must be taken to get unambiguous results. In this limit, 
one is faced with the problem that more and more cigcnvahics of C^^* get closer and closer 
to as y — 7> +00, which makes it difficult to tell them apart from ni and may jeopardize 
condensation by leading to a concentration onto a larger space than the eigenspace associated 
with K,i. For a homogeneous field^, one expects that issue to be all the more acute as the 
density of states at large wavelengths, close to the ground state (here, the condensate), is 
large. This will be the case at low space dimensionality d. Such a dimensional effect is well- 
known in traditional Bose-Einstein condensation of an ideal Bose gas which needs c? > 3 to 
exist. 

It is useful to summarize our main results. In this paper, we derive the criterion for 
the concentration of the Gaussian field (f (real or complex) to turn into a Bose-Einstein 
condensation as the thermodynamic limit is taken. The Gaussian field is characterized by 
its covariance function and we have the intensity u — WfWl/V as the control parameter. 
We show that in any dimension d (including, in particular, d = 1), there exists a class 
of covariance functions for which the system exhibits a Bose-Einstein condensation as one 
increases u through a critical value Uc. In this case, we investigate the precise "shape" of the 
condensate, i.e. the probability distribution of the part of u contained in the condensate. 
We show that depending on the parameters characterizing the covariance function and the 
spatial dimension d, there can be two distinct types of condensate: the first one is Gaussian 
distributed with normal fiuctuations scaling as 1 / VV ( "normal" condensate) , and the second 
one is not Gaussian distributed ("anomalous" condensate). A schematic phase diagram is 

^ i.e. with correlation function C{x,y) = C{x — y). Here and in the following, we use the word "homoge- 
neous" with the meaning of statistically invariant by translation in real or Fourier space, depending on 
the context (this will be specified explicitely in the text in case of ambiguity). 
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presented in Fig. 1. A detailed analysis of the structure of the anomalous condensate is 
performed for a one-dimensional system {d — 1). Extending this one-dimensional analysis to 
exactly the point of transition between normal and anomalous condensations, we find that 
the condensate at the transition point is still Gaussian distributed but with anomalously 
large fluctuations scaling as \/ln.{L)/L, where L is the system length (in this sense it can 
be termed "anomalous" too). 

We will see later that when expressed in terms of the Fourier components of the field, our 
problem is structurally very similar to the one-dimensional mass transport model introduced 
in [6], in particular the condensation properties in this model [7, 8]. There are however a 
couple of important differences. In the mass transport model, condensation happens in 
real space [7, 8], whereas in our Gaussian-field model, condensation happens in Fourier 
space. Secondly, unlike in the mass transport model where condensation in real space occurs 
homogeneously, i.e. the condensate can form at any point, here condensation is heterogeneous 
in Fourier space: for the particular class of fields we consider, it forms only at the A; = mode. 
This is more in line with the traditional Bose-Einstein condensation (see e.g. the review [9] 
for several other examples of homogeneous vs. heterogeneous condensation). Nevertheless, 
much of our present analysis can be performed along the same line as in [7, 8] in all the 
regimes, with and without condensation. Indeed, "normal" and "anomalous" condensates 
are also found in the homogeneous mass transport model [7, 8] , although the precise nature of 
the condensate found here differs from the one in the homogeneous mass transport model. It 
is worth mentioning that condensation has also been studied for zero-range types of processes 
on imhomogeneous networks [10], which is also structurally somewhat similar to our model. 

The outline of the paper is as follows. In Section II we specify the class of ip we consider 
and we give some necessary definitions. Section III deals with the conditions for occurence of 
Bose-Einstein condensation. The conditional spectral density of (f, knowing ii=||(^|||/V^, is 
given in Section IV for all the regimes (with and without condensation) . Finally, the struc- 
ture of the condensate is studied in Section V for a one-dimensional system in the different 
condensation regimes and extension to higher space dimensionality is briefly discussed. 



4 



II. MODEL AND DEFINITIONS 



Let A be a bounded subset of R''. In the following we take for A a d-dimensional torus 
of length L and volume V — L'^. The results in the thermodynamic limit are expected 
to be independent of the boundary conditions and imposing periodic boundary conditions 
by considering a torus makes the calculations simpler without loss of generality. Consider 
a homogeneous Gaussian random field 93 in A with zero mean and correlation function 
C(p(p*{x, y) = C{x,y) — C{x — y), supplemented with C^^ = if is complex and normalized 
such that C(0) = 1. Since A is a torus, C{x) is also periodic with period L. For every n e Z^, 
define 

Ci(n) = / C(x)e-'^'^"-^/^d'^x. (1) 

J A 

As a correlation function, C{x) is positive definite, hence CL{n) > 0. Moreover, C{x) is 
assumed to be such that (i) CL{n) < Cl{0) for every n 7^ and every L large enough, (ii) 
for every k eM.^ the limit 

C{k)^ lim CLin), (2) 

Z/,n— >-+oo 

27rn/L->fe 

exists and C{k) < (7(0) for every k 0, and (iii) C{k) = C{\k\) with C{k) ~ C{0) (1 - a|fc|") 
for small k (with a, a > 0). 

For every n G Z*^, we define the Fourier coefficients of ip as 

ipL{n) = ^ / ¥?(a;)e-2'™-"/^d'^a;. (3) 

W J A 



If (fi is complex, the (fiLin) are independent and the joint probability density function (pdf) 
of their real and imaginary parts is given by 

\Mn)r 



p[{Re(fiLin), ImifiLin)}] = TT J' exp 



Cr(n) 



(4) 



The interested reader will find the derivation of Expressions (4) and (5) in Appendix A. If 

(f is real, the <^L(n) are independent inside one given half of Z"* but not from one to the 
other complementary halves, as they are linked by Hermitian symmetry (pL[—n) — (fiin)*. 
It follows that a real (p is determined by half as many degrees of freedom as a complex </?: 
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its Fourier coefficients in only one fialf of Z*^, witfi joint pdf 

1 



p[{(pL{0),'R.e^L{n),lui(pL{n)}] = 



X 



V27rCi(0) 

n J- 



exp 



7rCi(n) 



exp 



2Cl(0) 



(5) 



"1/2 



wliere Z^^2 is a given fialf of excluding the point n — 0. Note that </7l(0) is necessarily 
real by Hermitian symmetry. Now, it is possible to recast ||v'||2 as a sum of independent 
random variables, which will prove very useful in the following. By Parseval's theorem one 
has 



I<^ll2 



(6) 



If Lp is complex, its Fourier coefficients are independent and (6) is already a sum of indepen- 
dent random variables. If (f is real, one rewrites (6) as 

Ml = M0r + 2 Yl [^^Mn)r + 2 i^^Mn)?, (7) 
which is also a sum of independent random variables. Thus, one can always write 



Ml 



(8) 



where the m„ are independent random variables given by quadratic forms of Re</?i(n) and 
Im (fL{n) according to either (6) (complex </?) or (7) (real </?). Namely, in the former case 
'iTT'n — |¥'l(^)P; whcrcas in the latter case one takes mo — </?l(0)^, = 2[Re (fiL{n)]'^, and 
rri-n = 2[lm(pL{n)]'^ for n e Since the m„ are independent, their joint pdf is a product 
measure straightforwardly obtained from either (4) (complex ip) or (5) (real (p). One finds 



P[{mn}] = Yl Pn{mn), 



with 



Pn[mn) 



4m„>0} 



TT 



lie 



exp 



mr 



(9) 



(10) 



where e — \ (resp. e — when is complex (resp. real)^. Replacing then ||(/?||2 with (8) in 
the expression p(t/) = (^dl^^Hi — TJ)) of the probability density of Hv^Hl, where the average 



^ A word of caution for the reader: the notation e should not be confused with the one of a small number. 
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(•) is over the product measure (9), one obtains 

p{u) = / ^u^-xi"^on ^^"("^") (11) 

J{mn>0} \ .^^a J riel^d 

At this point, it may be interesting to shortly digress and point out in more detail the 
similarities and differences between our model and the mass transport model studied in 
Rcfs. [6-8]. In the latter, there is positive mass variable m„ at each site n G Z of a one- 
dimensional lattice with periodic boundary conditions. The microscopic dynamics consists 
in transferring (or chipping) a random portion of mass from site n to site (n + 1). The 
amount of mass, m, to be transferred from a site with mass m to its neighbor is chosen from 
a prescribed distribution K {m\m) which is called the chipping kernel. The chipping kernel 
is homogeneous, i.e. does not depend on the site. The dynamics conserves the total mass 
M = '^n'lTin. For a class of chipping kernels [6], the system reaches a stationary state in 
the long time limit where the joint pdf of masses becomes time independent with the simple 
form [6] 

(M) = / 5 ( M - V m, I TT f(mn) drun, (13) 

■h'nn>0} \ i J n 

is the normalizing partition function, and N is the number of lattice sites. The weight 
function f{m), the same for each site n, is non-negative and depends on the chipping kernel 
K {mlm). By choosing this kernel appropriately, one can generate a whole class of weight 
functions /(m). The choice /(m) ~ m~'^ for large m with 7 > 2 leads to condensation 
whereby a finite fraction of the total mass M condenses onto a single site of the lattice when 
the density p = M/N exceeds a critical value [7, 8]. 

The product measure structure in the two problems is similar and p{U) given by (11) 
is just the analogue of the partition function Z^{M) given by (13). However, there is an 
important difference. In our case, the counterpart of the weight function in (11), Pni'iTT^n), 
depends explicitly on n (through its dependence on the covariance function CL{n)), whereas 
in the mass transport model the weight function finin) is identical for all sites n [6-8]. 
In this sense, our model is heterogeneous in contrast to the homogeneous mass transport 
model. 



P 

where 
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After this digression, we return to the expression (11) oi p{U) and express it in a 
more tractable form. The m„-integrals in (11) are most conveniently performed by us- 
ing the Laplace representation of the delta function or, equivalently, by writing p{U) as the 
Bromwich integral of its Laplace transform p{X) with 

r+oo 

p(A) = / p{U)exp{-XU)dU 
Jo 



-A 5^ m J \ 



= (^exp l-xymn]) (14) 

where we have used (8) and performed the average over the product measure (9), [recall that 
e = 1 (resp. e = 2) when ip is complex (resp. real)]. Then, by inverse Laplace transform, 

p{U)= ! e''U-\Y.^^^,H-^+eXCUn)\^^ (15) 
Jc 2i7r 

the integral being evaluated along a vertical line in the complex plane upward and to the 
right of all the singularities of the integrand. This contour of integration (or any of its 

continuous deformation within the analyticity domain of the integrand) defines a Bromwich 
contour generically denoted by C Making the shift s = A + 1/£Cl(0), one gets 

p{U)^e~^^) f e'^--e^n&dHi-CL{n)/CUO)+esCUn)]^_ (16) 

Jc 2i7r 
In the thermodynamic limit, it is more convenient to work with the intensive variable u — 
W^Wl/V instead of the extensive one U = \\ip\\l- Defining 

CL{n 



S{s, u,V) ^ su- 



cm 



+ £sCi(n) 



(17) 



and using (16), one finds that the probability density of u is given by the integral represen- 
tation, 

piu) = VpiU = uV) = Ve~^) I Qys{s,u,v)^_ (18) 

Jc 2z7r 

The leading asymptotic behavior of p{u) in the large V limit can then be obtained from a 
steepest descent analysis of (18) in which S{s, u, V) is approximated by 

su — I(s) (V ^ +00, s fixed), 

S{s,U,V)-^{ r -, ^ ^' (19) 

su- I{s) - esC{0)f{V)\ (s ^ 0, then!/ ^ +oo). 



where 



In 



1 



Cjk) 
C{0) 



+ esC{k) 



(20) 



and f{V) > 0, in the second hne of Eq. (19), is a non trivial correction resulting from 
the discreteness of the infrared modes around k — 0. It turns out that in the sequence of 
limits s — > then V +oo, these modes cannot be described properly as a continuum, 
and a careful analysis of the sum in (17) is required, leading to the correction fiV). In 
one dimension, f{L) can be obtained explicitely from the Euler-Maclaurin formula at lowest 
order applied to (17). One finds /(L) - with a = limfe^o In [^(A;) - C{Q)\/]n\k\. For 
higher dimensions, algebraic growth of f{V) is also expected as V +oo. The specific 
expression of f{V) will not be needed in the following. 

It is convenient to express the function I(s) in (20) in the following form 

d'^k 



I{s) = 7(0) + ^ 1 In (^1 + 



g{k)J (27r)'^' 



with 



9ik) 



1 



C{k) C(0). 

For later purposes we also define a critical intensity, Uc < +oo, by 



Ur 



/'(O) 



1 



(21) 



(22) 



(23) 



g{k) (27r)<^- 

The existence of the integral on the right-hand side of (23) depends on the small k behavior 
of g{k). For the class of ip we consider, one has 



£(7(0) 



(/c^O), 



(24) 



from which it follows that Uc < +00 (resp. Uc — +00) \i a < d (resp. a> d). 

The expression of p{u) depends on which limit in Eq. (19) must be used in Eq. (18). The 
right choice is determined by the contribution of the /c = mode to the sum (17) in the 
thermodynamic limit : if ti < Mc, it is negligible and the proper asymptotic expression for 
S{s, u, V) is the first line of (19) ; iiu > Uc, it becomes significant and the proper expression 
for S{s,u,V) is the second line of (19). (See the beginning of Sees. Ill A and IIIB, the 
discussion in Sec. VB, and Appendix B for further details). It is possible to write the 
corresponding two different expressions of p{u) into a single concise formula by introducing 
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a "tag", (T = or 1, indicating which hmit in Eq. (19) is considered. More precisely, injecting 
either hne of (19) into (18) one gets 

p{u) ^v\eC{0)f{V)y^\~^) [ s-<^Ae^[-«-^(*)]^, (25) 
L J 2nr 

where taking a = (resp. a — 1) gives the expression of p{u) obtained from to the first 
(resp. second) hne of (19). The criterion for a condensation of (f can then be found by 
analyzing the function 

5(ti) = - lim ^lnp{u). (26) 

A condensation transition of (p at some u = Ucond is indicated by the fact that 5(m) is not 
analytic at u — Ucond- Replacing p{u) with (25), one gets 

S(«) = _^- hm lln /s-/-e^[— ^W] — , (27) 
which is the expression of S(m) we will use in the next section. 



III. CONDITIONS FOR BOSE-EINSTEIN CONDENSATION 
A. u <Uc 

For fixed u < Uc there always is one saddle point s{u) > 0, solution tou — I'{s{u)), which 
determines the leading exponential behavior of the Bromwich integral on the right-hand side 
of (27). In this case S(m) reduces to 

E{u)^-^^+I{s{u))-us{u). (28) 

Note that s{u) being fixed, S{s,u,V) is given by the first Eq. (19), which corresponds to 
(T = Oin Eqs. (25) and (27). 

By analyticity of I{s) in C+ = {s e C : Res > 0} and the inverse function theorem 
for complex functions, s{u) is analytic in the interval < u < Uc- It follows that S(?i) is 
also analytic in the same interval and no condensation transition of (/? is to be expected 
for < w < Uc- Since Uc — +oo if o; > ci, one can already conclude that Bose-Einstein 
condensation of (p needs a < d to exist. 

We now take a closer look at the behavior of S(ii) ior a < d and u close to Uc < +oo. 
There are three different cases to consider: a < d/2, d/2 < a < d, and a — d/2. 
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1. a < d/2 and u u^. 



Write A{uy — —I"{s{u)) > 0. Prom the small k behavior of g{k), Eq. (24), it can be 
checked that A(0) < +oo if a < d/2. When u gets close to Uc from below, s{u) gets small 
and the saddle point equation can be expanded as 

u = I'{s{u)) ~ /'(O) + l"{0)s{u) 
= Uc- A{Oy s{u), 



yielding 

and, to the same accuracy. 



s{u) 



Ur — U 



I{s{u)) — us{u) 



m+i'io)siu)+-i"io)siur 



= /(O) + K - u)s{u) - -A{Of s{uf 



[UcS{u) + {u- Uc)s{u)] 



/(o) + 



2A(0)2 ■ 

Injecting (29) into the right-hand side of (28), one gets 

S(K) = /(0) + ^+^"'=-") 



£(7(0) 2A(0)2 ■ 



(29) 



(30) 



2. d/2 < a < d and u ^ 



In this case A(0) = +00 (because g{k) ^ is not infrared integrable) and a more careful 
analysis is needed. Take s > real and write (21) as 



I{s) = 7(0) +UcS+^ j 

gd/a r 

= 7(0) + UcS + — j 



In 1 + 



In 1 + 



g{k)) g{k)\ {2nY 
s \ s 



(31) 



d'^q 



gis^/'^q) J ^(s^/"g)J {2nY' 
where we have made the change of variable k = s^^°'q. Writing the last term at the lowest 
order in s, one gets 

I{s) = 7(0) + UcS - Js'^l'' + o(|s|'^/"), (32) 

with 



J 



£C(0) 



a\q\ 



- In 1 + 



a\q\ 



dU 
(27r)° 



(33) 



11 



Note that the inequahties d/2 < a < d ensure J to be both infrared and ultraviolet conver- 
gent. Since I{s) is analytic in C+, the expansion (32) holds not only for small s in R+ but 
also for small s in C"^. 

One can then follow the same line as in Sec. IIIAl with I[s) now given by (32). For 
u < Uc close to Uci the saddle point is solution to 



u = I'isiu)) c^Uc- —s{uY'''-\ 
a 



which gives 



s{u) 



a{uc — u) 
dJ 



and, to the same accuracy, 

I{s{u)) - us{u) ~ [7(0) + Ucs{u) - Js(ii)'^/°] - [ucs{u) + {u- Uc)s{u)\ 
= 7(0) + {uc - u)s{u) - Jsiuf'' 

Injecting (34) into the right-hand side of (28), one gets 

^/ \ T/^\ ( 0L\ (a\-^, ,_d_ 



(34) 



(35) 



3. a = d/2 and u u^. 



In this case the integral (33) is logarithmically ultraviolet divergent. To fix the problem 
we start like in the preceding section and split the /c-integral in the first line of (31) into the 
sum of one integral over |A:| > 5 and an other integral over |A:| < 6, where 5 > is arbitrarily 
small. Let A5(0) < -|-oo denote the regularization of A(0) = -|-cxd with an infrared cutoff at 
\k\ —5. Take s > real and make the change of variable k — s^/°'q in the low \k\ integral. 
One obtains 



I{s) ~ 7(0) + u^s - 

£ J\a 



MO? 2 

2 ' 
eC{0) 



(36) 



'|g|<(5/s2/d 



a\q 



I d/2 



- In 1 + 



i\q\d/^ j 



d'^q 
(2^' 



where the g-integral is now dominated by its ultraviolet behavior. Writing (36) at second 
order in s, one gets 



7(s) ~ 7(0) + UcS + Ks^ ln(s) - [As{Of + dK \u{5) + 0(1)] 
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2 ' 



(37) 



with 

^ - ^2^0^' ^^^^ 
where Sd is the unit sphere surface area in a d-dimensional space. Now, to deal with the 

logarithmic divergence of As{Oy as S tends to zero, we introduce a fixed r] > S small enough 

so that we can replace g{k) with its expansion near A; = for < rj, and we write 

MO) 



e g{ky (27r)<^ 

eC{QfSd p dk ^1 f 1 d'^k 



where 



{2nYa^ Js k ej^,\^^g{ky{2ny 



(27r)<^a2 ^'^ £y|,|>,^(A;)2(27r)<^' 
does not depend on S. It can easily be checked that dT/drj — too (as it should be, because 
A5(0)^ cannot depend on an arbitrary quantity like rj). For large |A;|, g{k)~'^ behaves like 
C{kY which is ultraviolet integrable, and T is finite. Using the expression (38) for K, we 
see that the logarithmic divergence of A5(0)^ compensate Kd\D.{5) exactly, yielding 

lim [A5(0)^ + dK ln(5)] = T < +oo, (39) 

and one gets 

7(s) =7(0)+'UcS + irs2ln(s) + o(|s^ln(s)|). (40) 

Since I{s) is analytic in C+, (40) holds not only for small s in IR+ but also for small s in C"*". 
For /(s) given by (40) and u < Uc close to Uc, the saddle point is solution to 

u — I'{s{u)) ~ lie + 2Ks{u) In s{u), 

which gives 

{u,-u)/2K 

s{u) 



\n{uc-u)/2Ky 
and, to the same accuracy, 

I{s{u)) — us{u) ~ [/(O) + Ucs{u) + K s{u)^ In s{u)] — [ucs{u) + {u — Uc)s{u)] 
= 1(0) + {uc - u)s{u) + (it) ^ Ins (it) 

c 1(0) + J- ~ . (41) 

4:K \ln{uc-u)/2K\ ^ ^ 
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Injecting (41) into the right-hand side of (28), one obtains 

H(„).,(0, + -4i^ + i,^-(jf|l^. (42) 

B. a < d and u> Uc 

In this regime, the contribution of the A; = mode becomes significant and the proper 
expression for S{s, u, V) is the second hne of (19), which corresponds to cr = 1 in (27). The 
size of the s-domain contributing significantly to the Bromwich integral on the right-hand 
side of (27) is then found to scale hke ~ 1/V{u — Uc) as V ^ +oo. For these s one has 
I{s) = /(O) + UcS + o{l/V) and 

/ „-l/e V[su-I{s)]_^ ^ -Vm f -l/e V{u-uc)s^ _ ^"^'^^"^ (A-X) 

from which it follows that (27) reduces to 

E{u) = /(O) + (44) 

C. Conditions for a condensation of 

If q; < d, it can be seen immediately from (30), (35), (42), and (44) that H(m) is not 
analytic at = Uc- This non analyticity indicates a phase transition oi aX, u — Uc- We 
will see in the next section that this phase transition corresponds actually to the onset of a 
condensation of </?, as mentioned at the end of Sec. II. 

On the other hand, \ia> d one has Uc = +oo and the saddle point equation u = /'(s(m)) 
has always a solution. No non analytic behavior of S(-u), is to be expected in this regime. 

The conditions for a Bose- Einstein condensation of ip are thus : a < d and u > Uc- A 
schematic phase diagram in the a-u plane (for a given fixed dimension d) is presented in 
Fig. 1. 



IV. THE CONDITIONAL SPECTRAL DENSITY OF THE FIELD (GIVEN u) 

To specify the nature of the phase transition undergone by (f at u — Uc ii a < d, we can 
analyze the conditional spectral density of (f for a given value of u and see what happens 
when u increases past Uc- 
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non condensed (or "fluid") phase 
05 

aid 



FIG. 1: Schematic phase diagram in the a — u plane for a given dimension d, where a > 
characterizes the small k behavior of the Fourier transform of the covariance function: C{k) ~ 
1 — a|/c|" as /c — )• (with a > 0). There exists a critical line Uc{oi,d) (solid line) separating the non 
condensed, or "fluid", phase (u < Uc{a,d), white region) and the condensed phase (u > Uc{o(,d), 
shaded region). The critical intensity Uc{a,d) diverges as a ^ d from below (vertical dotted line). 
For u > Uc{oi,d), there are two different types of condensates (separated by a vertical dashed 
line): a Gaussian distributed normal condensate for < a < d/2 and a non-Gaussian anomalous 
condensate for d/2 < a < d. 



Since u = Xlnez"^"/^ l^^^ (8)]' natural to define m„ = ntn/V as the contribution 
of the mode k = 27m /L to the intensity u. From the product measure (9) one can easily 
write down the conditional joint pdf of the random variables {u^}, conditioned to have a 
fixed sum u = Yln&'^ri- One finds 



P[{^n]\u] = ^^6[u-Y}\ \{Pr 



(45) 



where p(u) is given by (18) and, according to (10) and the change of variable m„ = rUnlV , 



Pn{Un) 



-{u„>0} 



7tV 



1/e 



1-1/e 
Un 



exp 



VUr 



(46) 



TT \eCLin) 

with e = 1, 2 corresponding respectively to the complex and real cases. The conditional 
marginal density of a particular Un is then obtained from (45) by integrating out all the 
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other variables Um^n- One gets 

P{Un\u) = , (47) 

p[u) 

where p±{u — Un) is the probabihty density of the incomplete sum Ylrnii^n '^"^ — u — Un (with 
given u and Un). Since the Um are independent random variables, p±_{u — Un) has the same 
functional form as p{u) in Eq. (18), except that S{s,u, V) is now replaced with 

Sn{s, u,V)^su ^7 In 



eV ^ 



(48) 



namely 

Jc 2z7r 

Using the conditional marginal distribution p{un\u) given by (47), we can define the condi- 
tional expectation 

E{Un\u) = / UnP{Un\u) dUn = —^-T / UnPn{Un)Pl.{u - Un) dUn- (50) 

Jo PW Jo 

It is easy to check that X]n^('"™l'") ~ should be. Physically, K{un\u) represents the 

average fraction of the total intensity u carried by the n-th mode. 

The conditional spectral density, Cu{k) — Cu{\k\), is subsequently defined in the ther- 
modynamic limit by requiring Cu{k) d'^k/{2T^)'^ to represent the average fraction of the total 
intensity u carried by all the modes that lie within an infinitesimal spectral volume d'^k 
around k. According to this definition, one has 

which is obviously normahzed to / Cu{k) d'^'k/ {27r)^ — u, as a, straightforward consequence 

of J2n'^i'^n\u) = U. 

Before proceeding with the determination of Cu{k) in the different regimes found in 
Sec. Ill, we first put E(m„|m) into a form more suitable to an asympotic estimate in the 
thermodynamic limit. From the expressions (18), (46), (49), and the Laplace transform 



lo " " e[{V/e){CM-^-CL{0)-^ + es)Y+^/^' 

(where we have written r(l -|- 1/e) — 7r^~^/^/£, valid for e — 1, 2), one finds after some 
straightforward algebra, 

E(u \u)-( f cMey^(--y) ds\ / f ^\ 

E{un\u) -[J^^_ c,{n)/Cm+esC,{n) 2in ) ^ A 2in J " ^^^^ 
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A. OL < d and u < Uc, or a > d: no condensation 



In this regime, the proper expression for S{s, u, V) in the large V hmit is the first hne 
of (19) and there always is one saddle point s{u) > 0, solution to u — I'{s{u)). This saddle 
point determines the Bromwich integrals on the right-hand side of (52) which reduces to 

hm y E(un\u)^ [ . (53) 

and, according to the definition (51), 

^ l-C{k)/C{0)+es{u)C{ky ^ ' 

Replacing I{s) with its expression (21) in the saddle point equation u = I'{s{u)) and using 
(54), one obtains J Cu{k) d'^k/ {27r)'^ = u, which means that (54) gives the complete dis- 
tribution of u among all the Fourier modes. It is clear from this result that, in the limit 
V — > -|-oo, the contribution of each given mode to u is infinitesimal and no condensate is 
present in this regime. 

Note that for u close to Uc, (or u large enough, if Uc — +oo), s{u) is small and (54) can 
be approximated by 

Cuik)^jJ^^ i\kr»esiu)/a), ^^^^ 
C'«(A:) ^ {\k\-<ss{u)/a), 
which depends on u in the vicinity oi k — only. 



B. cx < d and u > Uc'. condensation 

As mentioned in Sec. IIIB, the proper expression for S{s,u,V) in this regime is the 
second hne of (19). It follows that the size of the contributing s on the right-hand side of 

(52) goes to zero like ~ 1/V{u — Uc) as V" — )■ +oo, and for A; 7^ one has 



|n|>|fc|L/27r 

which yields 



lim > ¥.{un\u) = I ~ , , ,, 

^-+-,.,>1^.,. iM>|.|l-C(g)/C(0)(27r)^' 



C'(g) d^q 



(56) 



C (k) = ^^^'^ (57) 

' i-c{k)/c{oy ^ ' 

Replacing g{k) with its expression (22) in the definition (23) of Uc and using (57), one obtains 
Jkj^o Cu{k) d'^k/ {27ry — Uc < u. The normalization J Cu{k) d'^k/{2T^Y — u demands then 
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that the missing part u — Uc accumulate into the A; = mode. Let us check it exphcitely. 
For A; = the numerator on the right-hand side of (52) is proportional to 

(58) 

with £ = 1,2. And since the Bromwich integral in the denominator is proportional to (43) 
with the same coefficient of proportionality, one obtains 

lim Wuq\u) — u — Uc, (59) 

as expected. Combining Equations (56) and (59), one finds that the spectral density (57) 
can be prolongated to A; = 0, in the distributional sense, as 

Cu{k) = {2nf{u - Uc) 5{k) + J3^^^y^l{M0}- (60) 

The condensation into the A; = mode manifests itself in the presence of the 5-function on 
the right-hand side of (60). The sudden appearance of a condensate as u increases past Uc 
is a clear manifestation of the non analyticity of S(u) [or, equivalently, of p{u)] aX, u — Uc- 

V. THE STRUCTURE OF THE CONDENSATE {a < d AND u > tic) 

In this section, we investigate the structure of the condensate in more detail in terms of 
the asymptotic behavior as ^ -|-oo of the conditional marginal density of the zero mode 
uo (given u), 

p{uq\u) = — , 61 

p(u) 

for u > Uc- Using (18), (46), and (49) with n = in (61), one obtains 

p(^uo\u) = , ( f 5V^e^[^(^'"'^)-"°^l —)( f e^^(^'«'^) — V\ (62) 

Since u > Uc, the proper asymptotic expression for S{s, u, V) in the denominator of (62) is 
the second line of (19), and from (43) one gets 

p{uo\u) ^ V[eC{0)f{V)Y/^ (^)'~'^'"'''^°^ 

X f ,l/e^VlS(s,u,V)-uos]^_ (V ^ +Oo) (63) 

Jc 2nT 
The asymptotic behavior of the Bromwich integral on the right-hand side of (63) depends 
on the considered regime non trivially, as we will now see. 
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A. a< d/2 



Since Uq is expected to be close to u — Uc, the integral in (63) is expected to be determined 
by the contribution of the small s region. The actual size of this region determines which 
asymptotic expression (19) is to be chosen. Leaving this choice open for the moment and 
using (63), one obtains 



u — u, 



Uo 



1-1/e 



X 



f ^(l-a)/.gVW«-«o)-/(.)]^^ (V^+oo) (64) 

Jc 2z7r 



where a — (resp. cr = 1) if the first (resp. second) line of (19) is the right choice. Since 
A(0) < +00 if q; < d/2, one can write, for small s, 

I{s) = 7(0) + UcS - Ja(0)^s^ + o(|s|^). (65) 

Injecting (65) into the right-hand side of (64) and writing 

'= Aioy ^^^^ 

one gets 

(oi 7/ \ ^ y(uq—u-\-Uc)^ 



The value of a can now be determined by comparing, for the contributing s, the speed at 
which sCl{0) and 1 — Ci(|n| = 1)/Cl(0) in (17) tend to zero as 1/ — >■ +00. For instance, a 
rapidly decreasing sCl(O) going to zero faster than 1 — C^dnl = 1)/Cl{0) is equivalent to 
letting s ^ first, and then V +00. In this case, the asymptotic behavior of S{s,u, V) 
is given by the second line of (19) and one must take a — 1. In the opposite case where 
1 — CL{\n\ — 1)/Cl(0) goes to zero faster than sCl(O), one must take a — 0. So, we need to 
estimate the typical size of the contributing s. Because of the exponential in the first line 
of (67), the distribution of Uq is peaked around Uq = u — Uc with a width ~ l/\/V. Because 
of the exponential in the second line of (67), the contibuting t are in a small region around 
t — with a width ~ 1/VV. Thus, according to (66), the small s contributing to (64) for a 
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typical uq behave like 1/VV — l/L'^/^ as y — >■ +00, and since a < d/2, one finds 



lim 



C{0) 



lim sL" = 0, 



-Ci(|n| = l)/Ci(0) 

yielding cr = 1, as explained above. The integral in (67) reduces then to a simple Gaussian 
integral. Using {u — Uc)/uq — )■ 1 as y — > +00 (again, because the distribution is peaked 
around uq = u — Uc with a width ~ \/\/V), one obtains 



p{uo\u) 




_ y("0-" + "c) 



(y ^ +00) 



(68) 



27r A(0) 

In this regime, the part of u which is in the condensate is Gaussian distributed around u — u. 
with a width ~ 1/ \/V . We will call such a condensate a "normal condensate" . 



B. d/2 < a < d : statement of the problem and solution for d = 1 

The most natural way of proceeding to the case d/2 < a < d would be to start again 
from Expression (64) and determine the asymptotic behavior of the Bromwich integral by 
replacing I{s) with the small s expansion (32), instead of (65). Unfortunately, this line turns 
out to be a dead end. The size of the s contributing to (64) for a typical uq is found to 
behave hke l/V'^/'^ = for large V. It follows that (9(1) as V ^ +00 for the 

contributing s, and no definite value can be attributed to a. 

As we will see in the following, this unexpected problem is to be attributed to the fact 
that the asymptotics (19) are not correct if L^s — )■ 0(1) or L^s too slowly (e.g. 
logarithmically) as \^ ^ +00. A more careful, and finer, analysis of the asymptotic behavior 
of (17) is therefore needed to determine the one of p{uq\u) from the expression (63). 

We have been able to carry out such an analysis successfully for d = 1, taking advantage 
of the Euler-Maclaurin summation formula to control (17) and obtain p{uq\u) as a Laplace 
integral which can be computed numerically. The result is the existence of an "anomalous 
condensate" ii d/2 < a < d and u> Uc- 



1. Rewriting S{s,u,V) using the Euler-Maclaurin summation formula 
Prom the Euler-Maclaurin summation formula at lowest order, 

+00 ffn'l r+oo r+00 

n=l ^ -^0 
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where Pi{t) — t — n — l/2iorn<t<n + l, one gets 

+00 



El" 





f + OO 




In 




-00 


+ 2 


es — 





cm 



+ esCL{t) 



dt 



cm 



+00 



Plit)C'Lit) 



1 - 



cut) 

cm 



+ esCL{t) 



dt. (69) 



Making the change of variable t — A;L/27r, one obtains 



+00 



1 - 



cm) 
cm 





r+00 


+ esCm) 


= L In 




J —00 



cm) 



esCm) 



dk 
2^ 



+ 2 



es 



cm 



cm 



r-too 
Jo 1 - 



CLiO) 



+ esCm) 



-dk, (70) 



where we have written CmLf^n) = and C£(A;L/27r) = {2n/L)C'm)- Let R{s,L) 

denote the second term on the right hand side of (70). When s — > 0, i?(s, L) is logarithmically 
infrared divergent. To deal with this divergence, we rewrite R{s, L) as 

[Pi(i) + 1/2]C'1(A:) 



R{s,L) =ln esCm 



+ 2 



es — — 



cm 



CUk) 
Cl{0) 



+ esCm) 



-dk. 



(71) 



For = 1, it follows from (17), (70), and (21) that 



S(s,u,L)^su-I(s)~ ^^^'^^ (L^+00). 

eL 



(72) 



Our objective is now to determine the behavior of (72) as s — > with L°'s — > 0(1), and 
use the result in the Bromwich integral on the right-hand side of (63) to get the asymptotic 
behavior oi p{uo\u) as L — )■ -\-oo for d = 1, 1/2 < a < 1, and u > Uc. 



2. Determination of p{uq\u) for d = 1, 1/2 < a < 1, and u > Uc 

After some careful calculations detailed in Appendix B, one finds that in the limits 
L — >■ +00, s — >■ 0, and < 0(1), R{s, L) behaves as 

°r(l + l/a)r(l-l/a)f?5t(0)£^!V'°. (73) 
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with 

i2S r / f \ i^ 

(74) 



n=l ^ ' 



where C(-) is the Riemann zeta function. The second hne of (19) corresponds to the first term 
on the right-hand side of (73) only, with f(y) — L°'/a. Injecting (73) into (72), replacing 
I{s) with the small s expansion (32), and introducing the rescaled variables 

la yzvr / 

(2^^,_^ ^ 

£C(0) ^ ' 



one obtains p{uq\u) for uq close to u — Uc as 

»+oo 

£(7(0) 

Note that by the identity 



p{u,\u) r e--*-f^(-)|. (L ^ +OC) (75) 

£(7(0) J_oo 







1 , / 1 

In 1 



c^i = -ar ( 1+ - I r (i - - 

a I \ a 



the term in the second line of (73) exactly cancels out the term —Js^/"' from (32). The 
simple asymptotics given by the second line of (19) misses this cancellation. 

The ^-integral on the right-hand side of (75) can be computed numerically for fixed 
X. Then, varying X one obtains the kind of curves shown in Fig. 2. The non-Gaussian 
nature of the condensate in this regime is easily seen for, e.g., a = 0.7 and a = 0.9. More 
quantitatively, the skewness increases from 0.22 for a = 0.6, to 0.40 for a = 0.7, to 0.67 for 
a — 0.9. Similarly, the excess kurtosis increases from 0.11 for a — 0.6, to 0.39 for a — 0.7, to 
0.81 for a — 0.9. Note that, as a function of uo, the probability distribution p{uo\u) always 
tends to a delta function located ai uq — u — Uc + 0~ as L — > -|-oo. To see the structure of 
the condensate it is necessary to zoom in on the delta function, which is achieved by using 
the rescaled variable X instead of Uq, and rescaling p{uq\u) accordingly. The shifting of the 
peak of p{uq\u) to the left with increasing a, observed in Fig. 2, must be attributed to the 
fact that a — dis the limit for Bose-Einstein condensation (here d = 1). The closer to o; = o? 
one is, the larger L must be for the peak of p{uo\u) to eventually get close to uq — u — Uc 
by less than a given finite amount. The necessary L diverges in the limit a — > which is 
the precursor of the destruction of the condensate at a — d. 
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cL "{uq-u+Uc) 



FIG. 2: Resettled p{uq\u) vs. the rescaled variable X with c = (27r)"a/£C(0) (see the text for the 
definition of X), for d = 1 and a = 0.6 (solid line), a = 0.7 (dashed line), and a = 0.9 (dot-dashed 
line). 



In summary, at least for d = 1 where we are able to analyse the distribution p{uo\u) 
accurately, it turns out that the anomalous condensate in our case has a much more com- 
plex structure than the relatively simple form of the anomalous condensate found in the 
homogeneous mass transport model [7, 8]. 



3. Determination of p(uo\u) for d = 1, a = 1/2, and u > Uc 



The existence of an anomalous condensate for 1/2 < a < 1 raises the question of the 
structure of the condensate at the transition a = 1/2 between the normal and anomalous 
regimes. In this case, the behavior of R{s, L) in the limits L +00, s ^ 0, and L^/^s < 
0(1), is found to be given by 

R{s, L) - R{L^l^s) = In ( — ) + 2V^ ' ^ ' 



av27r 



+ 



eLKs"^ 



1 - 21n 



a\/27r 



with K = £C'(0)V(7ra2) [see Eq. (38)], and 



2 

n=l 



In ( 1 + — - — + — 
^^n J In 



(76) 



(77) 
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where 7 is the Euler-Mascheroni constant. The interested reader is referred to Appendix B 
for the details of the calculation. Injecting (76) into (72) with 



2 



I{s)c^I{0)+UcS + KsHn\^-^^j -(T+Kf-, (78) 
[see Eq. (37) with 0(1) made exphcit, and Eq. (39)], and introducing the rescaled variable 

ia \ 27r J 
one obtains p{uo\u) for uq close to u — Uc as 

p[uo\u) ~ — -1= / e at 

7rv2 A/, 



00 

2 



_ (mo~"+"c) 

1 e '""l 



/27r 

where 



(L ^ +00) (79) 



1/2 



a \ 7tL 

Like in the normal regime {a < 1/2), the part of u which is in the condensate is Gaussian 
distributed around u — Uc but with anomalously large fluctuations scaling as ^/\n(L)/L, 
instead of 1/ y/L for the normal condensate. 

It is interesting to notice that the size of the s contributing to (63) for a typical Uq in this 
regime is found to behave like l/^/L\nL for large L. It follows that L^/^s ~ l/^/\nL — )■ 
as L — )■ +00 for the contributing s, which suggests intuitively that the second line of (19) 
might be the right asymptotics to use. This is not the case. Like in the previous regime. 
Expression (79) results from non trivial cancellations between terms of (76) and (78) that 
are not accounted for by the simple asympotics (19). 



VI. SUMMARY AND PERSPECTIVES 



In this paper, we have investigated the concentration properties of a Gaussian random 
field ip in the thermodynamic limit, V +00 with fixed u — where V is the 

volume of the box containing the field. Considering a wide class of fields for which the 
spectral density behaves like C{k) ~ (7(0) (1 — a|A;|°) for small k (with a, o; > 0), we have 
found different regimes depending on the relative values of a and the space dimension d. 
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If q; < d and u > Uc, where Uc is given by Eq. (23), the concentration of (p discussed in 
[1] and [2] turns into a Bose-Einstein condensation onto the k — mode as the thermody- 
namic hmit is taken. If a < d/2, the condensate is normal in the sense that the part of u 
which is in the condensate is Gaussian distributed around u — Uc with normal fluctuations 
scaling as l/^/V. \i d/2 < a < d, the condensate is not Gaussian distributed ("anomalous" 
condensate) . A detailed semi-analytic analysis of the structure of the anomalous condensate 
has been performed ior d — 1. Extending this one-dimensional analysis to the transition 
between normal and anomalous condensations, we have found that the condensate at the 
transition is still Gaussian distributed but with anomalously large fluctuations scaling as 

\i a> d^ the density of states at large wavelengths gets too large as the thermodynamic 
limit is taken. The system cannot tell the infrared modes with > apart from the k — Q 
mode, which prevents Bose-Einstein condensation and leads instead to a concentration onto a 
larger function space than the only k — Q mode as is increased. This is the thermodynamic 
limit of the one-dimensional problem considered in [1], in which one typically has a — 2. 

For each regime, we have determined the conditional spectral density giving the average 
distribution of u among the different Fourier modes of the flcld. If a > d, or a < and 
u < Uc, the conditional spectral density is a smooth function of k given by Eq. (54). On the 
other hand, ii a < d and u > Uc, the conditional spectral density is given by Eq. (60) which 
clearly shows the condensation into the k — mode by the presence of the (5-function on its 
right-hand side. 

We have also illustrated the similarities and differences between the condensation prop- 
erties in the Gaussian-fleld model and those of the mass transport model in one dimension. 
While in both cases, there can be both normal and anomalous condensates depending on 
the parameter regimes, the precise nature of the condensate is different in the two cases. 
This difference can be traced back to the fact that while the condensation in the mass 
transport model is homogeneous (in real space) , it is heterogeneous (in Fourier space) in the 
Gaussian-fleld model. 

In conclusion we outline some possible generalizations of this work. As far as scalar flelds 
are concerned, further investigations would involve relaxing some of the assumptions we 
made about the correlation function. For instance, it would be interesting to investigate 
how our results are modifled when C{k) is not a single bump function of k. In the cases 
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where condensation occurs, the condensate is expected to have a much richer functional 
structure than a single Fourier mode, as it should live in the whole function space spanned 
by all the Fourier modes for which C{k) is maximum. Such a study would be of interest in 
e.g. laser-plasma interaction physics in the so-called "indirect-drive" scheme [11]. In scalar 
approximation, the crossing of several spatially incoherent laser beams leads to precisely this 
kind of Gaussian random field as a model for the laser electric field. An other interesting 
line of investigation would be to consider vector Gaussian random fields and quadratic forms 
more general than the simple L^-norm. For instance, in the context of turbulent dynamo 
[12], it would be of great interest to be able to determine the structure(s) of the random 
flow with a large average hclicity over a given scale. Assuming a Gaussian random flow, the 
question would then be to flnd out whether there is a concentration onto a smaller flow space 
when average helicity gets large and whether this concentration turns into a Bose-Einstein 
condensation when the considered scale gets large. Such studies will be the subject of a 
future work. 
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Appendix A: Joint probability distribution function of Re (pL{n) and Im<^i(n) 



If (p is complex with {<f{x)) — {(p{x)(p{y)) — and {(p{x)(p{y)*) — C{x — y), it follows 
from the definition (3) that the ^L{n) are complex Gaussian random variables with 



and 



{ipL{n)(pL{my 



{(fiL{n)) = {(fiL{n)(pL{m)) = 0, 
^ T? f [ M^MyT) e~^'''^'""''^d'^x(fy 

^ J A J A 

^ h [ [ C(a; - 7/) e-^^'^"-^^-^/)^^^, 

V J A Ua 



2iiT{'m—n)-y 



y 



(Al) 



(A2) 
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where we have used the periodicity of C {x) resulting from the periodic boundary conditions 
imposed on (/? and the definition (1) of CL{n). Writing (fiL{n) — Re(fL{n) + ilm(fiL{n) and 
Vhijn) = ReipL{m) + ilmipL^m) in (Al) and (A2), one obtains 



{Re(fiL{n)) = {Im (fiL{n)) = 0, 
{Rc ^L{'n)lnnpL{m)) = 0, 
{Re(pL{n)Re(pL{m)) = (Im(^L(n)Im(^i(m)) 



(A3) 



CL{n) 



from which it follows that the Re(pL{n) and the lmipL[n) are independent Gaussian random 
variables with zero mean and variance CL{n)/2. Their joint pdf is thus given by the product 
measure 

p[{ReipL{n), lmipL{n)}] = JJ A/'[o,Ci(n)/2](Re(^L(n)) x JJ A/'[o,Ci(n)/2](Im(^L(n)) 



= n 



1 



TTCLin) 



exp 



\Mn)\' 
CL{n) 



(A4) 



where A/[ci,c2](^)) with C2 > 0, is the Gaussian measure of mean ci and variance C2, 



\/27rc2 



exp 



[X - Ci, 

2co 



(A5) 



If (fi is real, its Fourier coefficients are linked by Hermitian symmetry, ipL{—n) — ipL{n)*, 
and (p is entirely determined by giving the (fLin) in only one half of Z''. Let 'Zf^^ be a given 
half of Z'^ excluding the point n — 0. All the statistical properties of cp are thus encoded 
in the joint pdf of (Pl{0) (which is necessarily real by Hermitian symmetry), Re(fiL{n), and 
Im ipiin) for n in 'Zy^. By ((/^(x)) = 0, {ip{x)ip{y)) = C{x — y) and the definition (3), the 
(piip) for n in {0} U Z^y2 ^'^^ complex Gaussian random variables with 



{^dn)) = 0, 



(A6) 



{(fL{n)(pL{m)) 



I f {v{x)ip{y)) e-'''<''-''+'^-y^d'^xd'^y 

^ J A J A 



J C{x-y)e-^'^''<''-y^(f. 

Cl{^ I ^ ^^^^^^ 

^ J A 



CLiO)5n,oSm,0, (A7) 
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{Sn-m — ^nfi^mfi because both n and m are in {0} U '^^2)^ 

= ^ f ^2Mm-n).y^ay ^ ^^^^^^^^ 
^ J A 



2in{m—n)-y 



y 



(A8) 



Writing (pL{n) — Re(fL{n)+ilm(fL{n) and (fiLiiTi) — Re</?L(m) + ilm</7i:,(m) in (A6), (A7), 
and ( A8) , one obtains 



(Re(/7i(n)) = (Im(/7i(n)) = 0, 

{Re(pL{n)lmipL{m)) = 

C (fi) 

{Re(pL{n)ReipL{m)) = ^ Sn,m, 



(A9) 



2-5. 



n,0 



for n and m in {0} U follows that <^l(0), the Re(/9i,(n), and the lm(pL{n) (for n in 

'^^2) independent Gaussian random variables with zero mean, ((/^^(O)^) = Cl{0), and 
((Re(^L(n))^) = {{Im ipL{n)y) = CL{n)/2. Their joint pdf is therefore given by the product 
measure 

p[{(Pl{0), Re(pLin), Im9?i(n)}] = 



\/27rCL(0) 



exp 



V^l(O)^ 



X 



n 



1 



7rCi(n) 



exp 



\VL{n)\' 
CL{n) 



(AlO) 



/2 



Appendix B: Calculation of R{s, L) 

We want to estimate the expression (71) of R{s, L) in the limits L — >■ +00, s — >■ 0, and 
< 0(1), with 1/2 <a<l. Since s tends to zero, on can already simphfy (71) as 

2 /•+°° [Pi{kL/27r) + 1/2]C^(A;) 



i?(s,L) =ln £sCl(0) 



(Bl) 



Cl{0) Jo 1 - CL{k)/CL{0) + esCL{k) 
For any given S > arbitrarily small, we spht the /c-integral on the right-hand side of (Bl) 
as the sum of an integral from to 5 plus an integral from 5 to +00. 
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The asymptotic behavior of the integral from 5 to +00 is easily obtained by noting that, 
in this domain, the derivative of ln[l — Ciik) / Cl{0) +esCL{k)] appearing in the integral is a 
well behaved function of k when s ^ and L — > +00. Thus, because of the wild oscillations 

oi Pi{kL/2n) + 1/2 around 1/2, one gets 



Cl(0) 



- In 



/■+00 
Js 1 - 



[Pi{kL/27T) + 1/2]C'M 



CL{k)/CL{0)+esCL{k) 



dk 



r+oo ^ 

Js 1 - 



C'{k)/C{0) 



C{k)/C{0) 



dk 



C{5) 
C{0) 



Inar. (L ^ +00, s ^ 0) 



(B2) 



In Eq. (B2), as well as in Eqs (B7) and (B9) below, the condition L"s < 0(1) supplementing 
the limits L — )■ +00 and s — )■ is always implicitely assumed. 

Since 5 is arbitrarily small, we can replace CL{k) with C{k) ~ C{0){1 — ak") in the 
integral from to S. Choosing S such that 5L/27r is an integer, one obtains 

[P,{kL/2n) + l/2]C',{k) ^;.^o, [Pi(A;L/27r) + l/2]afc"-^ ^^^ 



CLiO)Jo l-CL{k)/CL{0)+ssCLik) 



aL 



k^ 



<5L/27r-l 



TT Jo + 



dk — 2a 



n J 2: 



2n{n+l)/L j^a- 



+ rjs 



dk, 



(B3) 



n=0 -^STrn/L 

where we have used Pi{kL/27r) + 1/2 = kL/27r — n for 27rn/L < k < 27r(n + 1)/L and 
written rj — sCL{0)/a. We will now determine the asymptotic behavior of (B3), first for 
1/2 < a < 1, then for a = 1/2. 



a. 1/2 < a < 1 



The behavior of the first term in the second line of (B3) is given by 
aL r „ aSL 



aL r 

TT Jq /c" 



+ r]s 



dk 



aL 



TT 



S-T{1+ 



a 



TT 



r 1 



l-F 

1 
a 



f 1 1 


5" 


1, -, 1 + 




\ o; a 


rjs 




a 




— r)s 
a 



(B4) 
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where F{a, b, c, z) denotes the hypergeometric function 2-^1(0, b; c; z). Performing the inte- 
grals in the second term of the second hne of (B3), one gets 

SL/2TV-1 ^27r(n+l)/L 5L/27r-l 



2a n 

n=0 
<5L/27r-l 



+2 E 

n=0 

aSL 



n 



27rn/L + 

L°'r]s 



n 



n=0 



■ 27r(n + 1) , 2?™ 
m m —— 



In 1 + 



[27r(n + 1)]'* 



- In 1 + 



5L/2-K 



{2'KnY 

SL/2-K 



(B5) 



= ln(5-2Q;Vln — + — lnl + ^ -2Vln l + —\- 

n=l n=l \ \ / 

It is convenient to introduce the function 

+00 



n=l 



(B6) 



in terms of which one has the asymptotics 



(5L/27r 



n=l 



1-0 



(27r)" y ^ 2^ 1-a 



?7S (L +00, s 0), (B7) 



where we have used ES^''^"" ^ C(a) + (27r/5L)"[((5L/27r)/(l - a) + 1/2] (L ^ +00). 
Putting (B7), ln(l + ris/5") ~ 775/^" (s ^ 0), and 

(5L/27r 



2a ^ In 

n=l 

q;5L 



27rn a(5L, 27r , f5L\, 
ln^ + 2aln(^-j! 



In5 + aln^L 



TT L 

a5L 



TT 



TT 



(L ^- +00) 



(B8) 



in the place of the corresponding terms in the last line of (B5), one obtains 



SL/2n-l 2.(n+l)/L j^a-l 

2a y n dk ~ — 

±rf. J2nn/L k'^ + 7)8 TT 



n=0 

aln 6L + 20 



:l-o 



(27r)" 



1 — q; 
(L ^ +00, s ^ 0) 



(B9) 



Finally, it follows from (B2), (B4), and (B9) that the asymptotic behavior of (Bl) in the 
limits L +00, s ^ 0, and L"s < 0(1) is given by the expression (73). 
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b. a = 1/2 



The behavior of the first term in the second hne of (B3) is now given by 



dk 



r]s 



SL 
2^ 



1-F 1, 2, 3, 



27r 



5-2M^ln(^)-25VV' 



rjs 
(s^O) 



(BIO) 



and the Eqs. (B6) and (B7) are respectively replaced with 

+00 



m = c 



t 



111 

2 



n=l 



In 1 + 



t 



— + — 



(Bll) 



and 



5L/2-K 



^In 1 + 



n=l 



LV2 



TjS 



7rn 



rjs 



SL/2TT 



r]s 



27r 



271 



4:71 



5L/21T 



7 



(B12) 



n=l 



where we have used Ent/r^~^^^ ~ C(l/2) + {27i/5Lf''^{5L/7: + 1/2) (L ^ +oo) and 
Ent/r^"^ ~ 7 + ln((5L/27r) + 7r/(2(5L) (L ^ +cx)). As previously, In Eq. (B12), as well 
as in Eqs (B13) below, the condition L^^^s < 0(1) supplementing the limits L — > +oo and 
s — )■ is imphcitely assumed. Putting (B12), (B8) with a = 1/2, and ln(l + ~ 
■qs/S^/'^ — {rfsY / {25) {s — >■ 0) in the place of the corresponding terms in the last line of (B5) 
with q; = 1/2, one obtains 



5L/27r-l 

J2 n 



n=0 



dk 



27r(n+l)/L ^ 
27rn/L Vk + 7]S ^/k 



ln5L + 2^ 



TjS 



+ ^ [5 - 2S'/'ris - {vsf] + ^ ivsf In (^^) . (L ^ +oo, s ^ 0) (B13) 

Finally, using (B2), (BIO), (B13), and K = £(7(0)7(770^) = v'^/(s7r) [see Eq. (38)] in (Bl), 
one finds that for a = 1/2 the asymptotic behavior of R{s, L) in the limits L — >■ +oo, s -t- 0, 
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and L^/^s < is given by the expression (76). 
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